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Abstract. We study a natural implementation of Asymmetric Dark Matter in Twin Higgs 
models. The mirroring of the Standard Model strong sector suggests that a twin baryon 
with mass around 5 GeV is a natural Dark Matter candidate once a twin baryon number 
asymmetry comparable to the SM asymmetry is generated. We explore twin baryon Dark 
Matter in two different scenarios, one with minimal content in the twin sector and one with 
a complete copy of the SM, including a light twin photon. The essential requirements for 
successful thermal history are presented, and in doing so we address some of the cosmological 
issues common to many Twin Higgs models. The required interactions we introduce predict 
signatures at direct detection experiments and at the LHC. 
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1 Introduction 

There is overwhelming evidence for the presence of Dark Matter in the Universe, however its 
nature is yet to be unveiled. One particularly suggestive observation is that Ddm — 
raising the question of whether a common origin for the baryon and DM abundances is 
possible. One answer to this question is the Asymmetric Dark Matter (ADM) picture [1- 
6] (for recent reviews see [7-9]), which assumes the DM density to be determined by an 
asymmetry noM hr the same way the baryon number asymmetry ub sets the visible matter 
density. In general the abundances are related by 

_ "n-DM m-DM 
\Ib nB m-Tv 

where tudm is the DM mass and ttiat is the nucleon mass. If the two asymmetries are generated 
by the same mechanism or if one is responsible of the other then we expect udm ^ nB- It 
follows that mDM ~ 5m7v and the question to address is why this mass ratio is close to unity. 

The discovery of the Higgs boson and the lack of new physics signals at the LHC made 
Naturalness and the hierarchy problem a critical issue more than ever before. In search of 
solutions, renewed attention has been recently brought to Twin Higgs (TH) models [10-14] 
that solve the little hierarchy problem by introducing a copy of the SM. The mirroring of the 
SM Lagrangian is due to a Z2 symmetry. The scalar potential has an accidental global SU (4) 
symmetry which protects the Higgs potential from quadratic corrections. The Higgs is indeed 
a pseudo-Nambu-Goldstone boson and it is naturally light even after symmetry breaking, 
which happens at loop level. For Naturalness reasons the order parameter / of the SU{A) 
breaking is expected to be / ~ few v, with v = 246 GeV the SM Higgs vev. 

Twin Higgs is a natural environment in which to implement the ADM idea. It is straight¬ 
forward to draw a comparison between baryons and twin baryons and to consider the latter 
as DM candidates. Both are stable thanks to conservation of the SM (twin) baryon number 
B (B). Approximate Z 2 symmetry suggests udm ~ nB, for example by having only the 
combination B — B conserved above the scale /, and ttt-dm ~ with the second easily 
perturbed by the breaking of Z 2 below /. Motivated by this tantalizing observation we study 
the viability of twin baryons as ADM. 
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Related work [15-17] on both thermal and asymmetric DM has been recently carried out 
in the context of Fraternal TH models [18], in which only the third generation is mirrored. 
Here instead we focus on the conventional TH scenario with a copy of all generations. 


2 Twin sector properties 

Our starting point is a Twin Higgs model with an effective cutoff A r\j inf ~ 5 — 10 TeV, 
where again / ~ few v for Naturalness. We will suppose that there is a copy of each SM quark 
charged under twin SU{3) QCD and twin SU{2) weak gauge symmetries, whose couplings are 
expected to be close to the SM counterparts gi{K) ~ fl'i(A) (from now on all tilded symbols 
refer to twin sector quantities and fields). Hypercharge is not necessarily gauged, as recently 
shown in [18, 19], and eventually the twin photon is either massless or very heavy fn^ ~ A. 
The presence of leptons in the spectrum is not strictly necessary and we will comment on 
different possibilities in the next section. 

Naturalness constrains the twin top and bottom Yukawas, with yt and expected to 
fall respectively within 1% and 10% of the SM values, so from now on we fix yt = yt- While 
in general yt ~ yt is still expected from Z 2 symmetry, naturalness does not impose any bound 
on the first two generations. Hence large deviations can be possible, but we will not consider 
this in the following. Similar considerations apply to the leptons when present; either they 
will be assumed to have mass rht — f /v mi, or to be integrated out fhi 3> v. 

As for the interactions between the two sectors any TH model contains Higgs interactions 

CDyfhff-^^h^hff, ( 2 . 1 ) 

where h is the SM-like Higgs and its interactions with twin particles are those responsible for 
cancellation of quadratic divergences. At low energy they induce an effective operator 

■ ( 2 - 2 ) 

mhJ 

Higher dimensional operators will also be introduced in the following, with an expected effec¬ 
tive scale M ~ inf. In particular we have in mind recent realization of TH in the composite 
holographic Higgs framework [19-21], in which effective operators are generated by the strong 
sector and by integrating out heavy resonances. 

Finally the twin sector respects the same accidental global symmetries as the SM: twin 
baryon number, lepton number and charge. We will neglect possible small breakings generated 
by higher dimensional operators and consider stable each of the lightest particles carrying 
global charges. 


2.1 Twin Dark Matter 

Due to conservation of twin baryon number B, as already explained, twin proton p and 
eventually n are stable on cosmological time scales. If in the early Universe comparable 
asymmetries are generated then ub — and 

D r> TB ]\J 

^ _ 

D B mi\i 

where again mj\f {mi\f) is the (twin) nucleon mass and thus 

~ 5 ■ 


(2.3) 
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(2.4) 




As in the SM rriAr is determined by the QCD conhnement scale Aqcd, unless deviations of 
various order of magnitude are allowed in the Yukawas. Then satisfying Eq. (2.4) requires 

Aqcd — 5 Aqcd ■ (2.5) 


How likely is it to happen? The running of ag plays a crucial role. First the different masses 
of the twin quarks change the thresholds in the j3 function. It is also possible to directly 
introduce a small difference 6as = ds(A) — a*(A), for instance this can be due to threshold 
effects in composite TH realizations, in which heavy resonances are expected to have different 
representations under the twin and the SM gauge groups. 

We compute the conhnement scale of twin QCD by running dig down from dvr/, with the 
/? function computed for 6 havours of mass rhq = //\/2 yq. This leaves Sag, f and rescalings 
of all Yukawas (apart the yt) as variables. The results are presented in Fig. 1 which shows how 
Aqcd/Aqcd — 5 is a typical result once a 0(10%) splitting of the UV coupling constants is 
introduced. The effect can be explained by a simple approximate equation obtained by hxing 


Vi = Vi 

~ ( zV^'^Exp [— (— _ 

Aqcd~\v) ^[9 V «.( A ) ag{A))\ 


( 2 . 6 ) 


The hrst term accounts for the difference in the quark mass thresholds and its evident sup¬ 
pression explains the negligible sensitivity to the choice of / and y. On the other hand the 
initial condition is exponentially amplihed by the running. We also show the hne tuning (FT) 
contours given by two loop contribution to the Higgs mass [18] 


q„2 a2 

6ml ~ * {ug - dig) , (2.7) 

TT^ 

which is of the same order of the standard FT ~ (//t>)^. Then f /v < 5 is expected to avoid 
FT larger than few %. 

Having successfully accounted for the DM to baryon mass ratio, how can their number 
density be related? Following the ADM picture [ 6 ] some mechanism connecting the baryon 
asymmetries in the two sectors should be efficient in the early Universe. First consider an 
effective operator connecting SM and twin singlets, Z 2 invariant and preserving only B — B 


O 


BB — 


Amn 


( 2 . 8 ) 


where again M ~ dvr/ and all helds are right handed. The operator flavour structure, 
embedded in the coupling glpF, has been introduced to evade constraints on DM lifetime. A 
structure of the form gl'jfA oc yiyj..., as in partial compositeness, would be enough to avoid 
any bound and to ensure that the operator involving third generation quarks is in equilibrium 
at temperatures f RT < dvr/. Indeed for 5 ~ 1 the operator is efficient above temperatures 
T ^ M and, because it conserves B — B, '\t will convert any excess in B into B and viceversa 
so that it enforces hs ~ ns- 

A different possibility is the sharing of asymmetries through non-perturbative effects, 
for example leptogenesis followed by generation of both B and B. This solution depends 
crucially on the UV completion of the model and it is beyond the scope of the present work. 

We conclude, also anticipating the results in Section 3, that our DM candidate is a 
twin neutron h with mass rhn — 5 GeV. The h is typically stable and it is the main DM 
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Figure 1. The region of parameter space where 6 > ^qcd/^qcd > 4 as a function of / and 
Stts = 5s (A) — as(A). In blue the case in which yq = in orange the case yq = 2yq. The solid 
line corresponds to ^qcd/^qcd = 5 computed using the approximate function defined in Eq. (2.6). 
Dashed lines depict the fine tuning associated with two loops strong contributions to the Higgs mass, 
for FT values A = 5, 10, 20. 


component due to charge neutrality of the Universe if the B asymmetry is the only asymmetry 
generated in the twin sector. We will study its experimental signatures in Section 4. Let 
us remark that DM self interaction scales like {Aqqd/^qcd)^ resulting in a cross section 
Cn/hf-n — 0.25 cm^ g“^ which is below the current bound afi/mn < 0.5 cm^ g“^ [22, 23]. 

3 Thermal History 

A successful thermal history must address two main points: the absence of other, possibly 
overproduced, relics and the influence on BBN and CMB of additional relativistic degrees of 
freedom. 

Notice that we cannot rely on low reheating temperatures. The Higgs portal of Eq. (2.2) 
will always keep the two sectors in thermal equilibrium above temperature T ~ few GeV, 
while the DM (and baryon asymmetry) are expected to be generated at T > u. As such we 
discard reheating as a solution and from now on we always consider the two sectors to be in 
equilibrium up to at least T ~ few GeV. 

In the following we will present the thermal history, including a possible twin nucleosyn¬ 
thesis phase, for two scenarios differing in content of the twin sector: 

• Scenario A: a minimal scenario with only twin quarks. 

• Scenario B: a complete copy of the SM content, including a massless twin photon. 

3.1 Scenario A 

We start with a minimal scenario with only the necessary degrees of freedom, which are the 
twin quarks. We consider either the twin photon heavy or the twin hypercharge not gauged. 
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In the latter case there is no gange anomaly and no leptonic field is reqnired in the former 
all charged leptons and nentrinos are taken to be heavy enongh to be integrated ont. We will 
comment on this choice at the end of the section. 

This scenario does not inclnde any (nearly) massless particle and has no impact on the 
CMB. Only twin pions as possible symmetric relics are left to stndy. Their mass is given by 




'/ Aqcd 


V A, 


— 500 MeV 


QCD 


f 

700 GeV J 


(3.1) 


and similarly their decay constant scales as — ^QCD/^Qcnfn — 450 MeV. The are 
stable as they are the lightest charged particles. On the other hand is a singlet and can 
decay to SM states, with its width T being the crncial parameter. We assnme it to be indnced 
by higher dimensional operators. 

To stndy the cosmological history of the twin-pions we write the corresponding Boltz¬ 
mann eqnations. We assnme that operators responsible for decay keep the twin pions 
in thermal eqnilibrinm via the processes vr SM —>■ vr SM. Neglecting the twin baryons, the 
Boltzmann eqnations are 

no -h 3Hno = -{T){no - ho) - 2{av){nl - n±n±) , 
n± +3Hn± = —{(Tv){n±n± — n^) , (3-2) 

where no and n± are the nnmber densities of and respectively, and h represents the 
eqnilibrinm distribntion. Here, angle brackets represents thermal averages and a is the tttt 
scattering cross section 


cj(7r7r 


vrvr 


1 {s-mlf I dmj 

167r/4 s V s 


(3.3) 


Upon inspection of Eq. (3.2) it is clear that if T ^ 2H(T = m.^-) the decay is efficient in 
keeping fr in eqnilibrinm and no relic density will be prodnced. For T <C 2H{T = ttItt), 
eqnilibrinm is not maintained and the tt will become overabnndant, with each degree of 
freedom having a freeze-ont density 


Y = 0.278/5r*5(r = m^) . 


(3.4) 


After freeze-ont the term proportional to h in Eq. (3.2) drops to zero. The first term 
on the right hand side describes the decay of ffo, depleting their nnmber, while the sec¬ 
ond term replenishes it as long as an ;§> T. Above a certain temperatnre T the scat¬ 
tering will be efficient, n± will track no and the total nnmber density will decrease as 
hi{T) = hj(rjo)Exp[—r/(677(r))]. On the other hand below T the 7f=*= nnmber density 
will freeze-ont while all vf^ will eventnally decay. A rongh estimate of T can be compnted by 
imposing 


{av)sYExp[-r/{6H)] 

H 


T=T 


~ 1 


(3.5) 


^It must be pointed out that at least one extra twin-S'?/(2) doublet is necessary because of the global Witten 
anomaly [24], which requires an even number of doublets. For simplicity such additional field is supposed to 
have mass above the v, for example m = yffy^2) with y ~ 1, and be negligible in our discussion. 
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where s{T) is the entropy density. The hnal density of is then given by 


a 


~ 2 


m.^(so// 9 c) Exp -T/{m{f)) 


(3.6) 


An upper bonnd on the tt width can be estimated by plngging the solntion of Eq. (3.5) in 
Eq. (3.6). Imposing a conservative limit VL-^h? <1/10 ~ 0.012, we obtain 

f > 10“2® GeV . (3.7) 

Moreover, if the lifetime is too long its decays conld inject entropy dnring BBN. Thns 
we reqnire l/T < 1 s, corresponding to 

f > 6.6 X 10-25 GeV . (3.8) 


The bonnd of Eq. (3.8) antomatically satishes Eq. (3.7) and snggests that no signihcant 
relic density of is left and DM can be considered as solely composed of twin baryons. 
It is interesting to notice that snch a late decay conld be linked, and possibly solve, the 
cosmological lithinm problem [25]. 

Finally it is possible that for snfliciently long lived pions the Universe enters an early 
period of matter domination with snbseqnent reheating once the pions decay. This wonld 
dilnte both baryon and DM asymmetries. While irrelevant for present discnssion, it is an 
important point to address in any UV completion. 

If twin hypercharge is ganged the vr^ decays in SM via anomaly and 7-7 mixing [26] and 
the width is^ 


64vr3 



(3.9) 


where rh^ ~ TeV is the twin photon mass. However the condition T 2H is never satished 
when compared with the experimental constraint of e < IQ-^-lQ-^ [27, 28]. Alternatively, 
and in the case of hypercharge not ganged the decay can be obtained by introdncing 
higher dimensional operators. In particnlar, an axial cnrrent term (l/M^)( 57 ^ 75 ^)( 57 ^ 755 ) 
indnces a decay throngh mixing with the SM mesons. Notice that snch operators involve 
isospin singlets and so we inclnde isospin breaking effects in the form of the n-fj' mixing angle 
sin0 ~ 10-2. The width ronghly scales as 


T ~ 


rhl 


sin2 6<rx , 


(3.10) 


where the hrst factor acconnts for the mixing between the pion and an X meson, while the 
mass ratio acconnts for the phase space difference. To avoid additional isospin breaking 
snppressions from the SM sector we consider mixing with the rj' meson and so Eq. (3.8) 
translates into 

M < 2.5 TeV . (3.11) 

It is clear that for m,r < m^f = 958 MeV the reqnired scale wonld be too small for any 
reasonable model. The price to pay for heavy tt is in hne tnning as oc and thns we 
expect f /v > 10 (or drastic changes in the Ynkawas). 

^We are thankful to Tongyan Lin for pointing out an error in Eq. (3.9) of the original version of the paper 
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Another possibility is a 4-fermion interaction between twin quarks and SM leptons. The 
width is 

f (TfO ^ U) = sin2 e (^1 - , (3.12) 

where appears because of helicity suppression and decay to muons is the dominant channel. 
The bound from Eq. (3.8) gives 

/ ~ \ 

(jm) ■ 

which is less stringent than Eq. (3.11), but the presence of such operators is less easily 
justihable. 

What if twin leptons are added? Apparently the stability problem would be elim¬ 
inated if vr^ decay is allowed. However I annihilate only through weak interactions 

to V, and a modihed version of the Lee-Weinberg bound [29] then requires > few GeV 
to avoid overproduction. It would clearly reintroduce the stable problem. On the other 
hand, it is an interesting observation that a heavy i could be a WIMP candidate, a possibility 
recently explored in [15, 16]. The presence of light neutrinos does not signihcantly change 
the TT^ decay width, cause of the mass suppression from chirality flip, and bounds from CMB 
could be possibly constraining. We do not pursue this direction any further. 

Finally let us briefly comment on the possibility of twin nucleosynthesis, which could 
involve p and n once the temperature of the Universe is around the twin deuteron D binding 
energy. Notice that h is stable thanks to phase space. The fusion processes can proceed only 
with a pion in the hnal state with pn,pp,nh —)• At temperatures much above the 

binding energy the ratio of nucleons to nuclei is set by the high entropy of the Universe and 
the abundance of nuclei is negligible. In particular, this is valid at T < rh.,^ when the fusion 
is not efficient thanks to Boltzmann suppression. We conclude that DM is entirely composed 
of h and possibly p (if hypercharge is not gauged). 


3.2 Scenario B 


Now we consider a complete copy of the SM and again assume y* ~ y*. Various experimental 
constraints [27, 28] require a small photon-twin photon kinetic mixing e < 10“®, which could 
be possible if mixing arises at 4-loop or higher order. A second constraint comes from the 
number of effective neutrino species, the recent Planck measurements [30] give AN^ff = 
0.11 ± 0.23. The number of relativistic degrees of freedom is dehned 


i=bosons ^ ^ 2 =fermions ^ ^ 


(3.14) 


where gi accounts for the internal degrees of freedom. In general in the twin sector is com¬ 
parable to the one in the SM. Above the decoupling temperature the only difference is given 
by the larger mass of twin particles, which become non relativistic at a higher temperature. 
After decoupling the different temperature of the two sectors must be taken into account. 

If the twin neutrinos are light, the contribution of D and 7 is equal to the SM one, 
5 * = 3.36, rescaled to take into account the different temperatures of the two sectors 


ANeff ~ 


3.36 

0I5 



(3.15) 
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where the normalization factor corresponds to g'* = 0.45 of one SM neutrino species. By 
conservation of entropy the temperatures ratio is 


T 

T 


9*{Td) y^^ 


9* {Td) J 


(3.16) 


with Td the temperature of decoupling of the two sectors. If the decoupling happens at 


the freeze-out of Higgs portal interactions, then ~ 9*{Td) and AN^ff ~ 7, which is 

ruled out. The best possibility is to have the sectors decouple between the two QCD phase 
transitions when git{Td) ~ 61.75 while git{Td) ~ 10.75. Thus we obtain AN^ff ~ 0.7, which 
is however still in tension with Planck measurements. 

The most natural solution is to decouple the u by lifting their masses, for instance by a 
suitable choice of the twin seesaw scale. The effect is two-fold: increasing the entropy ratio 
at decoupling, ~ 5.5, and decreasing the total number of relativistic species as can be 

readily seen by a modihed version of Eq. (3.15) 



(3.17) 


which is compatible with present constraints. We conclude that the two sectors must decouple 
between the two phase transitions with 0.15 GeV ^ ^ 1 GeV. How can this be achieved? 

Additional interactions between the two sectors must be introduced and e should be involved. 
For a 4-fermion operator suppressed by a scale M the decoupling temperature is roughly 
given by the condition 



(3.18) 


Then the required scale is 7 TeV ^ M < 30 TeV. We stress that the above results are 
independent of /. The estimate only involves rhe, which is few order of magnitudes lower 
than Aqcd for — He, and thus our conclusions are not sensitive to /. 

If the neutrinos are integrated out, then each leptonic species is stable due to individual 
lepton number conservation and could in principle be a thermal relic. However, the U —)• 77 
annihilation channel is efficient enough to eliminate any symmetric component [31]. 

In this scenario vf^ as well as re are stable due to phase space, while ffo will decay to 
77 through the twin chiral anomaly. The relic density of is negligible, as shown in the 
previous section. 

What about nucleosynthesis? First consider the case in which the only asymmetry 
generated in the twin sector is in B. Then charge neutrality of the Universe imposes that 
only re are generated after the twin QGD phase transition. Twin neutrinos are heavy and 
so weak processes of the form e~^h -f-)- rep are suppressed. We are left with fusion processes 
nh —)• DTr~ and as in the previous section we conclude that DM is completely composed of 
re. 

On the other hand, suppose a L asymmetry is generated in the early Universe as well. 
Now p and e are also present, and by charge neutrality hp = fie- So the fusion process 
ph — 7 - Dy is open and twin atoms will be formed. The ratio of different nuclei is set by the 
initial leptonic asymmetry and we leave a detailed study of this case for future work. 
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4 Experimental signatures 


The experimental signatures predicted by the model fall in two main categories: DM direct 
detection and possible collider signals. Both of them are crucially dependent on the higher 
dimensional operators that we have introduced in Section 3. In general the presence of 4- 
fermion operators induce a DM-nucleon scattering cross section. For simplicity we consider 
an operator of the form 

■ (4.1) 

As previously explained we assume that DM is in the form of twin neutron with mass — 5 
GeV. The spin independent cross section on a nucleon N is 

where ^ is the reduced mass while bp = + b^, bn = bu + 26^ and similarly for twin baryons. 

The parameter space with M < 2> TeV is already probed by present experiments [32-34] as 
shown is shown in Fig. 2 for flavour universal couplings bjsf = bn = 3. 

Notice that the lowest values of M are demanded for decay in Scenario A, and only 
axial couplings are strictly necessary. In that case only spin dependent are induced with cross 
section of magnitude comparable to Eq. (4.2) which is below present bounds. 

On the other hand, Higgs portal interactions will always contribute with 


asi{n-N) 


ifN'mN)^{fnrhn)‘^ 2 
A £d ^ 

TTmtP 


~ 10 


-48 


cm 


-2 


(4.3) 


where /„ ~ /^r ~ 0.3 is the nucleon matrix element. The predicted cross section is below the 
neutrino floor. 

The presence of effective operators will give additional collider signatures on top of the 
usual TH signatures [13]. Twin quarks are pair produced and any twin particle will typically 
escape the detector. Mono-jet [35, 36] and mono-photon searches [37] can be safely evaded 
for M > 1 TeV. For a more tuned scenario in which the vr are heavy enough to decay inside 
the detector it could be possible to have signatures similar to emerging jets [38]. 


5 Conclusions 

We have shown that Twin Higgs models provide a natural framework for Asymmetric Dark 
Matter. They automatically contain the required global symmetries and particles, in partic¬ 
ular considering the parallel between baryons and twin baryons. 

Indeed, twin baryons with mass ~ 5 GeV can be naturally obtained by small Z 2 breaking 
effects which induce a higher twin QGD confining scale. It is also natural to expect the same 
order of magnitude for the baryon asymmetries in the two sectors. We have studied two 
scenarios differing in spectrum and interactions, and shown that twin neutrons are typical 
DM candidates. 

Regardless of the twin particle content, higher dimensional operators with effective scale 
M ~ 0(1 — 10) TeV are expected and often necessary for a successful thermal history. Their 
introduction provides experimental signatures within current or near future reach, both at 
Dark Matter direct detection experiments and at the LHG. 
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moM [GeV] 

Figure 2. DM direct detection parameter space. The shaded region corresponds to DM-nucleon spin 
independent cross sections mediated by the operator Oi for 6 > tudm > 4 GeV and M > 1 TeV. The 
lines correspond to the upper bounds from current DM direct detection experiments: CMDSlite [32] 
(dotted blue), SuperCDMS [34] (dashed purple), LUX [33] (solid red). 
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